arXiv: 1502.00200v 1 [cs.NI] 1 Feb 2015 


Base Station Switching Problem for Green Cellular 
Networks with Social Spider Algorithm 

James J.Q- Yu, Student Member, IEEE and Victor O.K. Li, Fellow, IEEE 
Department of Electrical and Electronic Engineering 
The University of Hong Kong 
Email: {jqyu, vli}@eee.hku.hk 


Abstract —With the recent explosion in mobile data, the energy 
consumption and carbon footprint of the mobile communications 
industry is rapidly increasing. It is critical to develop more 
energy-efficient systems in order to reduce the potential harmful 
effects to the environment. One potential strategy is to switch off 
some of the under-utilized base stations during off-peak hours. 
In this paper, we propose a binary Social Spider Algorithm 
to give guidelines for selecting base stations to switch off. In 
our implementation, we use a penalty function to formulate the 
problem and manage to bypass the large number of constraints 
in the original optimization problem. We adopt several randomly 
generated cellular networks for simulation and the results indi¬ 
cate that our algorithm can generate superior performance. 

Index Terms —Green cellular network, base station switching, 
social spider algorithm, evolutionary computation, swarm intel¬ 
ligence. 

I. Introduction 

With the growing concern on the global climate change, 
it becomes critical to develop energy-efficient systems in all 
industries. In particular, the information and communication 
technology (ICT) sector was estimated to contribute to more 
than 830 million tons of carbon dioxide emission in 2013, 
which is approximately 2% of the global carbon dioxide 
emission, and this number is expected to double by the year 
2020 m. With the fast growing mobile data communication 
demand, the wireless cellular networks are playing a more im¬ 
portant role in ICT sector than ever before. Thus a promising 
approach to reduce the global green house gas emission is 
to reduce the energy consumption of the cellular networks. 
It is also very important from the economical perspective of 
network operators as a significant portion of the operational 
expenditure is due to electricity consumption. 

Designing green cellular networks, especially green base 
stations, is a recent hot research topic. There are at least two 
mainstream approaches. With the development of smart grid 
technology, one approach is to utilize renewable energy for 
base station operation Q. Another approach is to manage the 
operation profile of the base stations. As the network operators 
need to deploy their base stations to support the peak mobile 
data traffic, it is inevitable that during a major portion of the 
day a large number of the base stations are under-utilized. This 
phenomenon can be observed from real mobile data usage 
profile 0. Unfortunately, the energy consumption of an idle 
base station is nearly the same as one under full load j4j. So 
an under-utilized base station must be switched off in order to 


save energy, and there naturally rises a problem: Which base 
stations should be switched off to save the most amount of 
energy while maintaining adequate service for end users? 

Such energy-efficient design of wireless cellular networks 
has attracted significant attention recently. However, most of 
the existing efforts focus on the formulation of the problem 
00. This paper focuses on designing an swarm intelligence 
algorithm based on the Social Spider Algorithm (SSA) ID- 
In this paper, we concentrate on the modifications on the 
original problem to make it solvable with metaheuristics, and 
the implementation of a metaheuristic. 

The rest of the paper is organized as follows. In Section 
[II] we briefly introduce the current progress on green cellular 
network design. In Section [III] we elaborate on the original op¬ 
timization problem and our modifications. Section |TV| presents 
the detailed implementation of our proposed algorithm. Sec¬ 
tion [V] introduces the random network generation protocol, 
and the simulation results are presented with analysis and 
discussion. Finally we conclude in Section [VI] with potential 
future research. 

II. Background 

Green cellular network design is a recent hot research topic, 
and much work has been carried out. Chiaraviglio et al. 
initially proposed the idea of dynamic base station operation 
based on the data traffic profile in 0. Later they extended their 
work and set up an analytical model to control the base station 
switching profile 0- Oh and Krishnamachari also proposed a 
similar model in 0, and devised a simple base station switch¬ 
ing policy. They further developed a complete system model 
for this problem and proposed a distributed algorithm to solve 
it 0. In this work, both the algorithm design and the practical 
implementation, including a distributed base station switching 
protocols were elaborated, and the preliminary simulation 
results indicate a satisfactory energy-saving performance. All 
the above work focus on base station operations in the same 
type of networks, i.e. macro cellular networks. There are also 
some researchers concentrating on the cooperative operation 
considering network sharing. Fehske et al. investigated the 
possibility of deploying small and low power base stations 
alongside conventional ones from a deployment perspective in 
0- Mars an and Meo evaluated the energy savings achieved 
with the energy-aware cooperative management of the cellular 
access networks of different operators over the same area in 


ED- in HD, Rost and Fettweis discussed the energy-efficient 
operation based on the cooperative transmission in multi-hop 
systems. Niu et al. introduced the concept of cell zooming 
for energy savings. In their implementation, the cell size is 
adaptively adjusted according to different control variables. 


III. Problem Formulation 


In this paper, we adopt the system model formulated in 
ED as it is a simple yet generalized wireless cellular network 
model. In this model, the system load is formulated as follows: 



tOM) 

r b (x,t) 


dx, 


(1) 


where p b (t) represents the system load of base station b at time 
t, A b represents b ’s coverage, 7 ( 0 ;, t ) represents the traffic load 
of user equipment x at t, and r b (x,t ) represents the service 
rate of x from b at t. This system load stands for the fraction of 
time needed to serve all the data transmission load in the base 
station’s coverage. For simplicity, we omit the time variable 
in the following notations. 

With this system model, we further formulate our Base 
Station Switching Problem (BSSP) as follows. We consider 
that the system load for a wireless cellular network remains 
constant in a short time interval. An optimization problem may 
be formulated as follows: 


min E a b 

bGB , (2) 

s.t. 0 < p b < pb, \/b € B 


where B is the collection of all available base stations, a b £ 
{ 0 , 1 } is the active indicator of &, and ~p^ is the maximum 
allowed system load for b. When a base station b is switched 
off, its traffic load will be handled by all its active neighboring 
base stations N b = {N bt i, N b g, ■ ■ ■ ,N btTl }, and we use p b ^i 
to denote the amount of transferred traffic load from b to the 
*-th neighbor base station in N b . 

From ([2]) we can see the fitness evaluation function is very 
simple, but the rigid constraints may potentially obstruct the 
generation of feasible solutions to the problem. In order to 
alleviate the effort in designing an algorithm that can easily 
generate solutions satisfying all the constraints, we transform 
the constraint in Q into a penalty function. The new fitness 
function is described as follows: 


min ^(a b + \B\ x P b ), (3) 

beB 

where P b is the penalty value for b. In order to define P b , we 
first analyze the possible scenarios a base station may come 
across during the switching operation: 

1) & is switched off, and |iV&| = 0. In this case, both b 
and all its neighboring base stations are switched off 
and the traffic load originally handled by b cannot be 
served anymore. P b is defined as 1 + p b for all such 
base stations. 

2) b is switched off, and \N b | > 0. In this case, b is switched 
off, but it has some active neighboring base stations. The 


traffic load of b is handed over to TV*, and the transferred 
load to N b j is p b ^i. P b is defined as 0 for all such base 
stations. 

3) b remains active, and a number of its neighboring base 
stations are switched off. We use S b to denote the 
collection of the neighboring inactive base stations. Then 
the system load for & is p b + Y^ s eS b Ps^-b- If this 
new system load is larger than ~p b , P b is defined as 
1 + p b + J2ses - Pb- Otherwise, P b = 0 for all 
such base stations. 

This penalty function is designed according to the design 
principles that are critical for a successful penalty function 
stated in Q 2 : 

1) The fitness values shall improve as the solutions ap¬ 
proach feasible regions of the search space. 

2) The unfeasible solutions are guaranteed to be assigned 
with fitness values inferior to the fitness value of the 
worst feasible solution. 

Thus, any solutions with fitness values larger than \B\ are 
unfeasible ones to the original optimization problem ([ 2 ]). 

It is noted that BSSP is NP-complete as it can be reduced 
from a vertex-covering problem, which has been shown to be 
NP-complete lfl3l . 

IV. Algorithm Design 

In this section, we will first briefly introduce SSA. Then the 
detailed implementation of our proposed methodology will be 
presented. 

A. Social Spider Algorithm 

SSA is a recently proposed general-purpose swarm intel¬ 
ligence algorithm. It mimics the foraging behavior of the 
social spiders to perform optimization task. In SSA, the 
search space of the optimization problem is formulated as 
a hyper-dimensional spider web, and each position on the 
web represents a feasible solution. Besides the solution space, 
the spider web also serves as the transmission media of the 
vibrations generated by the spiders. 

The spiders are the basic operating agents of SSA. Each 
spider holds a memory consisting of its current position on 
the web, the fitness value of its current position, and the 
vibration (position and intensity) it was following in the 
previous iteration. The first two pieces of information describe 
the characteristic of this spider, and the last one helps the 
algorithm guide the movement of this spider. 

Based on biological observations, spiders are found to 
be extremely sensitive to vibrations. They can accurately 
sense the strength and the direction of vibrations, and can 
even separate different vibrations propagated on the same 
web fl4l . SSA utilizes this characteristic of the spiders and 
established a vibration-based information-loss communication 
system among spiders. In SSA, a spider generates a new 
vibration whenever it moves to a different position from the 
previous one. The vibration will then propagate over the spider 
web and be sensed by others. By this means, the spiders 



share their personal information to form a collective social 
knowledge. 

The vibrations are defined by two properties in SSA, 
namely, the source position and the source intensity. When 
a spider moves to a new position, it generates a vibration at 
its current position. We define the intensity of a vibration in 
the range [0, -Too), and the value of the intensity is calculated 
as follows: 

l/{Cmax ~ f(s)) for maximization 
1 /(f( s ) — C m in ) for minimization 

where / is the intensity of the vibration at its source position, 
f(s) is the fitness value of the source position, C max is a 
confidently large constant selected such that all possible fitness 
values of the maximization problem are smaller than C max , 
and G m i n is a confidently small constant such that all possible 
fitness values of the minimization problem are larger than 
Grain,' 

([4]) provides the method of calculating the source intensity 
of the vibration. However, as a kind of energy, the vibration 
attenuates over time and distance during the process of prop¬ 
agation. This physical phenomenon is also taken into consid¬ 
eration in SSA from two aspects, i.e., time and distance. As 
a swarm intelligence algorithm, SSA performs the searching 
task in an iterative manner. The vibrations are attenuated with 
each iteration as follows: 

I(t + 1) = 1(f) x r 0 , (5) 

where I(t ) is the source vibration intensity at iteration t, 
and r a is a user-defined vibration attenuation parameter. This 
design can prevent the algorithm from pre-mature convergence 
as a non-decaying vibration on the web can potentially attracts 
all spiders to move continuously towards it, thus facilitating the 
exploitation searching behavior but obstructing the exploration 
of the whole search space. 

Another attenuation factor besides the time is the propaga¬ 
tion distance. In SSA, the vibration attenuation over distance 
is defined as follows: 

I{p) = I(s) x exp(^ ), (6) 

where s and p are the source position and receiving position 
of the vibration, respectively, J(s) is the vibration intensity at 
s, D(s,p) is the distance between s and p, and D max is the 
maximum possible distance between any two positions on the 
web. We usually employ the Manhattan distance to reduce the 
computational time. 

There are three phases in SSA, namely, initialization, it¬ 
eration, and final. In each run of SSA, we start with the 
initialization phase where the objective function, search space, 
and the optimization parameters of SSA are initialized. Then a 
random population of spiders are generated and placed on the 
spider web, i.e., search space. The positions of these spiders 
are randomly generated in the search space, and the initial 
vibration each spider is following has a position of the spider’s 


current position, and an intensity of zero. This finishes the 
initialization phase and SSA proceeds to the iteration phase. 

In the iteration phase, the algorithm performs the optimiza¬ 
tion task in an iterative manner. A number of iterations are 
executed. In each iteration, the algorithm first evaluates the 
fitness values of all spiders on the web and attenuates the 
vibrations in the previous iteration. Then the spiders generate 
vibrations using (|4]) at their current positions. The vibrations 
then propagates using 0 over the spider web. After the 
propagation, each spider will receive a number of different 
vibrations from all directions. Upon the receipt of these 
attenuated vibrations, each spider chooses one vibration with 
the largest intensity, i.e., the strongest vibration Vb es t- Vb es t is 
then compared with the stored vibration this spider followed 
in the previous iteration v prev . The one with a larger intensity 
is kept and saved as Vbest■ The algorithm then manipulate the 
position of the spiders as follows: 

P s (t +1) = P s (f) + (Pbest (i) - P s (t )) © (1 - R © R), (7) 

where P s (t) is the position of spider s at iteration t, and 
© denotes element-wise multiplication. Pbest is the vibration 
source position of the best vibration Vb es t- R is a vector of 
random numbers generated from zero to one uniformly, whose 
length is the number of dimensions of the problem, and 1 is 
a vector of ones of the same length as R. 

After the manipulation of spider positions, an artificial 
spider jump away process is introduced in order to prevent 
the algorithm from getting stuck in the local optimums. Each 
spider has a small possibility of being re-assigned with a 
new random position in the search space. And this ends one 
iteration. 

The iteration phase loops until the stopping criteria is 
matched. The algorithm then outputs the best solution with 
the fitness value. The above three phases constitute SSA and 
interested readers can refer to ED for more details. 

B. The Proposed Implementation 

SSA was initially proposed to solve global optimization 
problems, i.e., the solution space is continuous. So we make 
several modifications to adapt it to solve BSSP, which is a 
discrete problem. 

1) Encoding Scheme: In this work, we adopt the classical 
encoding scheme for this kind of problem. We use a vector of 
0 ’s and l’s to represent the off/on state of each available base 
station. 

2) Spider Following Scheme: In SSA, we use 0 to ma¬ 
nipulate the movement of spiders. However this equation is 
designed for continuous optimization. So we devise a new 
spider following scheme to replace 0 as follows: 

P s (t + 1) = P s (t) + round((P bes t(t) - P s (t)) © R), (8) 

where round () is the rounding function. 0 will first de¬ 
termine the dimensions at which Pbest{t ) and Ps(t) are 
different. Then all different dimensions have a probability of 
0.5 to change. 



3) Spider Jump Away Scheme: Besides the spider following 
scheme, another scheme that can change the position of the 
spiders is the spider jump away scheme. In SSA, the jump 
away scheme operates at the spider level, where each selected 
spider is assigned with a new random position. But in our 
implementation for solving BSSP, the jump away scheme op¬ 
erates on the dimension level. Right after the spider following 
scheme has manipulated the position of each spider, each 
dimension of a spider position has a probability of l/|f?| to 
change from zero to one or one to zero, i.e., 

P s {t + 1) = \step{R) - P s {t)\, (9) 

where step(r) = 1 if r < 1/|.B|, otherwise step(r) = 0. 

These are the modifications we made to adapt SSA to solve 
BSSR Algorithm [T] is the pseudo-code of SSA for BSSP 


Algorithm 1 Social Spider Algorithm for Base Sta- 

tion Switching Problem _ 

1 : Assign values to the parameters of SSA. 

2 : Create the population of spiders pop and assign memory 
for them. 

3: Initialize v prev for each spider. 

4: while stopping criteria not met do 
5: for each spider s in pop do 

6 : Evaluate the fitness value of s. 

1: Attenuate the intensity of v prev . 

S: Generate a vibration at the position of s. 

9: end for 

10 : for each spider s in pop do 

11 : Calculate the intensity of the vibrations V 

generated by other spiders. 

12 : Select the strongest vibration Vb es t from V. 

13: if The intensity of Vb es t is smaller than v prev then 

14: Stole V prev as Vbest - 

15: end if 

16: Perform the spider following move towards v be st - 

17: Perform the spider jump away modification on 

each spider. 

IB: Store Vbest us V prev . 

19: end for 

20 : end while 

21 : Output the best solution found. 


V. Experimental Setting and Simulation Results 

In order to evaluate the performance of saving energy using 
our proposed algorithm, we performed a series of simulations 
over some randomly generated instances. 

A. Testing Instance 

In our simulation, we consider the 3G network topologies 
consisting of 20, 40, and 60 base stations in the area of 10 x 10 
km 2 . The instances are generated as follows. We first randomly 
deploy the base stations to different positions in the area. Then 


we start from the first base station and determine the number 
of its neighboring base stations using 

\Neigh b \ = \N b + S b \ = Pois {A - 2) + 2, (10) 

where Pois( A) is the Poisson distribution and A is the 
expected number of neighboring base stations. Each base 
station may have a different number of neighborhoods, and 
this design guarantees that all base stations will have at least 
two neighborhoods. 

After determining the number of neighborhoods for each 
base stations, we establish the neighboring links starting from 
the first base station to the last one. In this step we first 
check the existing links connected to the current base station. 
If the number is already larger than \Neigh b \, we go on to 
the next base station. Otherwise we establish a neighboring 
link between the current base station and its nearest base 
station that is not linked yet. This substep will iterate until 
enough links are established. The pseudo-code for this random 
instance generation algorithm is presented in Algorithm [2] 


Algorithm 2 Random Instance Generation Algo¬ 
rithm_ 

1: for each base station b in B do 
2: Assign a random position in the area for b. 

3: end for 

4: for each base station b in B do 

5: Determine the number of V s neighboring base stations 

\Neigh b \. 

6 : end for 

7: for each base station b in B do 

B: Find the number of b’ s established neighboring links 

n. 

9: if n > \Neigh b \ then 

10 : continue 

11 : else 

12: Find the \Neigh b \ — n nearest base stations to b 

that not yet linked with b. 

13: Establish neighboring links between these base 

stations with b. 

14: end if 

15: end for 

16: Output the best solution found. 


Using this method we generated 9 instances using the 
instance parameter (|.B|,A) £ [(20,3), (40,4), (60,6)], and 3 
different instances are generated for each instance parameter 
pair. These instances are shown in Fig. [T| where the dots are 
the base stations and the base stations connected with lines are 
neighboring base stations. The first two digits in the names of 
the instances are B\ and A, respectively, and the last letter 
stands for different instances with the same parameters. 

B. Experimental Setting 

To determine which base stations to switch off, we need 
some feedback information from all the base stations. However 
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(a) Mean results of 20 base stations with A = 3 
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Fig. 1. Plot of Randomly Generated Topology Instances with Different 
Parameters 


these feedback information may reduce the overall perfor¬ 
mance and sometimes is invalid. So in our simulation we use 
an approximation method to reduce the reliability of these 
external and internal information. We make two assumptions 
on the random instances generated: 

1) The system load for each base station is equal to the 
overall system load: pb = p s ,Vb £ B. Here p s is the 
overall system load. 

2) The transferred traffic load from one base station to all 
its neighbor base stations are equal: pb^i = Pb/\B\ 

With these two assumptions we minimize the impact of all 
additional information and focus on the optimization perfor¬ 
mance of our proposed algorithm. 

We compare the performance of our proposed SSA for 
BSSP with the SWitching-on/off based Energy Saving algo¬ 
rithm (SWES) El whose system model was adopted in our 
problem formulation. In particular, the SWES( 0 0 ) algorithm 
is employed for performance comparison as no additional 
information (external or internal) is required for this algo¬ 
rithm, thus satisfying our assumptions. For SSA, we use 
the number of base stations as the size of the population, 
and the vibration attenuation rate is set to 0.9. As the spi¬ 
der jump away scheme is modified, we do not need the 
jump away rate as elaborated in 0. The maximum itera¬ 
tion number is set to 500, so for 20-base-station instances 
the fitness function is evaluated for 10 000 times, pi, is 
set to 0.6 for all base stations in all the instances and 
we test the performance of the compared algorithms with 
Ps £ [0.05,0.1,0.15,0.2,0.25,0.3,0.35,0.4,0.45,0.5,0.55] 
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Fig. 2. Comparison of Active Base Stations 



for each instance. 

C. Simulation Results 

Fig. [2] illustrates the simulation results of our generated ran¬ 
dom topology instances with SSA for BSSP and SWES( 0 ,o)- 
Each dot in the figure is the mean value of the number of active 















































































Fig. 3. Plot of Randomly Generated Topology Instances with Different A 

base stations in all three instances with the corresponding same 
set of parameters. From the simulation results we have the 
following observations: 

1) SSA for BSSP always outperforms SWES( 0 ,o) in all 
three set of simulations. Out of 33 data points, SSA 
performs better than SWES(q,o) in 29, and the remaining 
4 comparisons end up with draws. 

2) Both SSA for BSSP and SWES(o,o) are efficient in 
terms of saving energy. Assume the traffic profile is 
increasing and decreasing evenly during the day, SSA 
can save 42.5% of the total energy and SWESjo.o) can 
save 31.3%. These numbers are generated by averaging 
all data points indicated in Fig. [2] 

3) The advantage of SSA for BSSP over SWES(o : o) is more 
significant when the system load is relatively low. A 
potential reason is that when the system load is close 
to the maximum allowed system load, the number of 
base stations that can potentially be switched off is very 
small. This means that the number of feasible solutions 
to (f2]i is very small and a heuristic like SWES( 0 ,o) can 
also find a good solution to the problem. In the extreme 
case of p = 0.55, no base station can be switched off 
and the performance of the two compared algorithm 
naturally becomes the same. 

D. Impact of A on Energy Saving Performance 

Besides the comparison of SSA for BSSP and SWES(o,o)> 
we also analyze the impact of A on the performance of saving 
energy. In order to have a complete analysis on this issue, we 
generated 5 topologies with A € [3, 5,6, 7,8], \B\ =40 and 
conduct simulations on these instances as well as the “40.4.A” 
instance illustrated in Fig. |T| The topology instances to be 
considered are presented in Fig. [3] 

The simulation results on these new instances with SSA for 
BSSP is shown in Fig. [4] From the results we can observe 
that although A = 3 performs slightly worse than others when 
the system load is high and low, the overall performance of 



System load p s 

Fig. 4. Simulation Results of Instances with Different A 

these instances are barely distinguishable. This phenomenon 
suggests that increasing the neighboring links among the base 
stations would not significantly improve the overall energy 
saving performance when the system already has sufficient 
links. It can also be observed that the performance curves are 
very close to the ideal best performing curve p s /~pb- This again 
shows the outstanding performance our proposed algorithm. 

VI. Conclusion and Future Work 

In this paper we consider the possibility of saving energy 
in wireless cellular network by switching off the under¬ 
utilized base stations during non-peak hours. We transform 
the original optimization problem with constraints into an 
unconstrained problem with a penalty function to alleviate 
the design difficulty considering the constraints. Besides the 
new formulation, we also propose a new binary SSA to solve 
this optimization problem. SSA is a newly proposed general- 
purpose metaheuristic. It was originally devised to solve con¬ 
tinuous optimization problem. We make several modifications 
to the canonical SSA to adapt it to solve BSSP. We devise two 
optimization schemes in SSA and make modifications on the 
control flow of the algorithm. We employ a series of randomly 
generated network topologies to test the performance of our 
proposed SSA for BSSP. The simulation results indicate that 
SSA can almost always outperform the compared algorithm. 
We also analyze the impact of the network density on the 
energy saving performance. The results imply that the perfor¬ 
mance improvement caused by the increasing network density 
is limited as the performance is already very close to the ideal 
best case scenario. 

There are some potential future work. In this paper, we 
consider the base stations consume the same amount of energy 
despite different system loads and communication ranges. In 
reality, both factors will influence the energy consumed. This 
can be taken into consideration but will significantly increase 
the complexity of the optimization problem. Another extension 
is to consider the diversity of the base stations. In this work we 
only consider homogeneous base stations. It will be interesting 


































to study the impact of adopting different types of base stations 
on the energy-saving performance. 
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